Stability of equatorial circular geodesies in static axially symmetric spacetimes 
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A general study of the stability of equatorial circular orbits in static axially symmetric gravitating 
systems is presented. Important circular geodesies as the marginally stable orbit, the marginally 
bounded orbit and the photon orbit are analyzed. We found general expressions for the radius, 
specific energy, specific angular momentum and the radius of the marginally stable orbit, both for 
null and timelike circular geodesies. Solutions expressed in cylindrical coordinates, oblate spheroidal 
coordinates, and prolate spheroidal coordinates are considered. We show that all null circular orbits 
are unstable and that there are not marginally stable null geodesies, whereas that for timelike 
geodesies the orbits can be unbounded, bounded or circulars. 
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I. INTRODUCTION 

A major problem in the General Theory of Relativ- 
ity is obtaining exact solutions of the Einstein equations 
corresponding to the gravitational field of acceptable con- 
figurations of matter. Now then, as one of the most fun- 
damental characteristic of isolated systems in universe 
is axial symmetry, static or stationary axially symmet- 
ric exact solutions are of great astrophysical relevance. 
Accordingly, through the years, a great deal of work has 
been dedicated to the theoretical study of this kind of 
exact solutions [l|. 

On the other hand, once the solutions are found, an- 
other equally important problem, closely related with the 
above study, is the analysis of the motion of test particles 
in the gravitational field generated by such distributions 
of matter. Indeed, the study of the motion of test par- 
ticles provides valuable information about the structure 
and behavior of such gravitational fields. Furthermore, 
the study of orbits in the equatorial plane is of clear astro- 
physical relevance due to its relation with the dynamics 
of intergalactic stellar motion or the flow of particles in 
accretion disks around black holes. 

The motion of test particles in axially symmetric space- 
times has been studied by different authors through the 
years, both for satic as for stationary spacetimes and 
with different configurations of sources (see, for instance, 
references Q to [3])- Now, the purpose of the preset 
work is a general study of the stability of circular or- 
bits in the equatorial plane in different gravitating sys- 
tems formed by axially symmetric structures. In partic- 
ular, we will analyze some important circular geodesies 
as the marginally stable orbit, the marginally bounded 
orbit and the photon orbit. 

The paper is organized as follows. Section|TT]is devoted 
to derive the geodesic equations, the effective potential, 
and general expressions for the main characteristic of cir- 



cular orbits: the radius, specific energy, specific angular 
momentum and the radius of the marginally stable orbit, 
both for null and timelike geodesies. Then, in the follow- 
ing sections, we particularize these expressions for solu- 
tions written in cylindrical coordinates, oblate spheroidal 
coordinates and prolate spheroidal coordinates. 

Thus, in section ITOI for solutions expressed in cylindri- 
cal coordinates we concluded that all null circular orbits 
are unstable, as is illustrated by considering the Chazy- 
Curzon field. Then, in section ITVl we present the oblate 
spheroidal coordinates and some members of the family 
of Morgan-Morgan disks are analyzed. Later, in section 
IVl we consider prolate spheroidal coordinates and the 
range of stability of the Erez-Rosen solution is obtained. 
Finally, results are discussed in section IVll 



II. TEST PARTICLE MOTION 

The metric for a static axially symmetric spacetime 
can be written as the Weyl line element [l|, 



ds 1 



-e 2 ^dt 2 + e- 2 ^[p 2 dip 2 + e 2 ^(dp 2 + dz 2 )}, (I) 



where 7 and tp are functions of p and z only. The ranges 
of the coordinates (p, tp, z) are the usual for cylindrical 
coordinates and — 00 < t < 00. The Einstein vacuum 
equations reduce to the system of Weyl equations [H, HU 
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where ([2]) is the well-known Laplace equation in cylindri- 
cal coordinates with axial symmetry, which is the inte- 
grability condition of the overdetermined system ((3])-(j4|). 

The corresponding Lagrangian for this line element ([!]) 
is given by 



2C 



^12 



'i 2 + e~ 2 ^[p 2 tp 2 + e^(p 2 + z 2 )}, (5) 
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where the dot represents the derivative with respect to 
the affine parameter along the geodesic, A. Now, as the 
Lagrangian is independent of t and ip, 

-E = dC/di, £ = dC/dip, (6) 

are conserved quantities, where £ is the specific angular 

I 



momentum and E is the specific energy with respect to 
infinity. 



From the Lagrangian ([5]) we can derive the system of 
motion equations 



p + (p 2 - i 2 )( 7 ,p - t P ) + 2pi(7,, " i>,z) + e- 27 
z - (p 2 - i 2 )( 7 , z - ^,») + 2pz( 7:(9 - ^ iP ) + e- 27 V, 



P A 



£ 2 e 4 ^ 



(7) 
(8) 



which has a unique solution when conditions Xq = x a (Xo) 
and Uq = i a (Ao) are given, with x a — p,z. The initial 
condition for the velocity of the particle is obtained by 
substitution of © in ([5]). 

Now, if we confine the motion of the particle to the 
equatorial plane z = 0, from ([5} we obtain for the radial 
coordinate p the equation 



•2 -2-i 

p = e 1 



E 2 - e 2 V - 



(9) 



Now, for circular orbits we have that p = p c = 
constant and so p = p = 0. Accordingly, from expression 
(HIJ follows that 



E z 



V(p), 



(14) 



with V{p) given by (|12[) . Furthermore, the minimums of 
V{p) correspond to stable circular orbits, whereas that 
the maximums of V{p) correspond to unstable circular 
orbits. So, by computing the derivative of V(p) we obtain 
the equation for their critical values, which can be written 



with e = 1 for timelike geodesies and e = for null 
geodesies. The orbit of the particle in the equatorial 
plane can be obtained by solving together the above equa- 
tion and the equation 



£e^L 



0, 



(10) 



that follows from For radial motion we have that 
tp = (po = constant, so that £ — in the above equations. 

The behavior of the trajectories in the equatorial plane 
is determined by the equation ([5]), that can be conve- 
niently expressed as 



p 2 e 2 ^ 



= 2V 



= 2V 



= — . (11) 

2 v ' 



so that we can define an effective potential through 



V(p) 



= e 2 ^ 



e 2 + V* 
P 



(12) 



which only depends on p and the metric function tp. On 
the other hand, in order that the metric ([1]) be asymptot- 
ically flat, the functions "0 and 7 must vanish at infinity. 
So, we can obtain the general condition 



lim V 

p— >oo 



(13) 



^ 2 e 2 ' / '(2p^-l) + p 3 e 2 ^ = (15) 

and, for the case of null circular orbits (e = 0), as 

2p^ p -1 = 0. (16) 

So, the radius of the timelike and null circular orbits 
are given, respectively, by the roots of the two previous 
equations. 

The specific angular momentum £ for massive particles 
in circular orbits can be obtained from equation (1151) and 
is given by 



p 3 ^, P e' 



-2ip 



1 " 2pt, 



(17) 



with the condition < pip^ p < 1/2. So, we can see that 
the radius of the circular orbits it depends on £. Now, 
by replacing (|17|) into equation (fT4|) . we obtain the other 
constant of motion, E, for a particle moving in a circular 
trajectory 



E* = 



5^(1 -W,p) 
(1-W,p) ^ 



(18) 



where, again, < pip tP < 1/2. 

The stability condition for circular orbits is given by 
V"(p c ) > 0. So, for massless particles the stability con- 
dition reduces to 



for all the effective potentials of the form (fT2")) . 



ijj. p + pip.p P > 0, 



(19) 



3 



whereas that for massive particles with specific angular 
momentum given by (|17p , the stability condition is given 
by 



pV, P p + 3V> + 2pV 2 p (2pV',p - 3) > 0, 

with < pip ;P < 1/2. 

Now, one can show that the expressions 



V"(p) = 



and 



dp 







(20) 



(21) 



(22) 



are equivalent. Accordingly, the radius of the marginally 
stable circular orbit can be obtained through the two si- 
multaneous equations V'(p) = and V"(p) = or by 
means of the equation dt 2 /dp — 0, provided that there 
exist two critical points of the effective potential, one 
of them corresponding to the stable circular orbit and 
the other one to the unstable circular orbit. Thus, the 
minimum value of the specific angular momentum as a 
function of the radius of the circular orbit (|17l) , it repre- 
sents the last circular orbit, which is well-known as the 
marginally stable circular orbit. For null geodesies the 
expression is 



ip,p+pip, PP = 0, 
and for timelike geodesic is 



(23) 



p V,p P + 3^ + 2 (0 </> 2 p (2p^ p - 3) = 0, (24) 

where < pip, p < 1/2 again. 

Finally, we can also find an expression for the angular 
velocity, 



dip cp 
~dt ~ 7' 



(25) 



wherein i = Ee 2 ^ and ip = £e 2 ^p 2 . For timelike 
geodesic we obtain 



(26) 



where p c are the roots of equation (|15p . whereas for null 
geodesies we have 



,2 V 



(27) 



III. SOLUTIONS IN CYLINDRICAL 
COORDINATES 

In spherical coordinates (r,9), the asymptotically flat 
solutions of the equations system ((2J - (j4]) are [l[ 



7fe 



r n+l ' 



n=0 
k 

E 

Lm=0 



CiC m {l + l)(m+l) 
(l + m + 2)r l+m + 2 



(28) 



(29) 



where P n = P n (cos 9) are the usual Legendre polynomials 
and the C n are constants. Now, in the equatorial plane 
z = so that we have 



ti'h 



Ik 



2^ ° 2 » pln+l 



(30) 



n=0 
k 

E 



C 2 iC 2m {2l + l){2m+l) 
(2l + 2m + 2)p( l + m )/ 2 + 1 



l,m=0 

X (P2lP2m - P 2 /+l^Wl)- (31) 

where (p, z) are the usual cylindrical coordinates, with 

p = r sin 9, z = r cos 9. 

We can obtain the radius of a circular null orbit by 
solving the equation (|16l) . that in these coordinates re- 
duces to 



E 

n=0 



2C 2n P 2 „(0)(2n+l) 



n 2n+l 



1. 



(32) 



So, if C 2rl P 2 n(0) > 0, we obtain 



2fc+l \ " 
n=0 



2n n 

p =o, 



(33) 



a polynomial in p off odd order that has, at least, a real 
root. Moreover, since there is only one change in sign, 
there is a positive root. Accordingly, we can conclude 
that there exist null circular orbits. On the other hand, it 
is easy to see that the stability condition it is not satisfied, 
since 



pipk, PP + V'fe,, 



E 



C 2n P 2n (0)(2n + l) 2 

p 2 n+2 \ > 



where p c are the solutions of expression (|16j) . The above 
equations only depend on the metric function tp, in such a 
way that the potential 7 it is not needed for a qualitative 
analysis of the particle trajectories in Weyl spacetimes. 
However, the function 7 is necessary for solving the dif- 
ferential equations of motion of the particle. 



and so all the circular orbits are unstable. Finally, we 
can ask for the existence of a marginally stable circular 



orbit, which must satisfy pip. 



pp 



0, but wc find 



that there are not positive roots as the corresponding 
polynomial has not any change of sign. Therefore, there 
is not any marginally stable orbit. 
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Now, in order to illustrate the above considerations, 
we take the simplest case of the family ([25)1 , the Chazy- 
Curzon solution frzllill. 



Wo 



rn 
r 



m 2 sin 2 i 



7o 



2r 2 



(35) 



which can be obtained taking k = and Co = m > 
in (|28p. As we can see, although the metric function ipo 
is spherically symmetric, the full solution (|35|) is not. In 
the equatorial plane, the metric functions reduce to 



m 

Wo = j 

P 



7o 



2p 2 



(36) 



So, according to (|3"2"j) . for this solution the radius of the 
unstable circular orbit is p — 2m. Then, with p = 2m, 
the values of the specific energy and the angular velocity 
of the particle are 



cjtv = (2me) 



E 



N 



£uj n , 



(37) 



where £ is an arbitrary constant. 

In Fig. [JJwe show the effective potential for lightlike 
geodesies in the Chazy-Curzon solution. As we can see, 
there is a maximum at p = 2m with a value of (2me)~ 2 £ 2 . 
Trajectories can be described using the horizontal lines 
(V = E 2 ) to different values of the quantities m and £. 
When < E\ < (2me)~ 2 £ 2 the motion corresponds to 
a particle with specific energy E\ coming from infinity 
until reach the turning point B and then going back to 
infinity. There is also a potential well, when the motion 
it is confined to < p/m < A, and a not allowed region, 
for A < p/m < B. On the other hand, when the specific 
energy it is greatest than (2me)~ 2 £ 2 as in the horizontal 
line E2, there are no turning points and the particle it 
moves only in one direction. Now, although we do not 
consider here solutions with other values of k in ([3H|) . it 
can be shown that in general all the effective potentials 
behave as depicted in Fig. [TJ whenever the condition 
C2n-P2n(0) > be assumed in this family. 

On the other hand, for timelike geodesies the specific 
energy and specific angular momentum of a particle in a 
circular orbit are, respectively, 



£ 2 = 



a 2m/p 



p — 2m 
where p > 2m, so that 



E 2 = e~ 2m/p , 



2 

0J T 



-Am/ p 



(p — 2m)p 2 



(38) 



(39) 



The radius of the marginally stable circular orbit can be 
obtained by solving the equation (j2"4"]l and so we can find 
the corresponding specific angular momentum by replac- 
ing this radius in equation (|17j). For the Chazy-Curzon 
field we have 



p = 5.23m, 
E = 0.826, 



= 3.52m, 
, 0.00525 



77 J" 



(40) 
(41) 



for the marginally stable circular orbit. 

The graphics of the effective potential (fT2")l for timclikc 
geodesies in the Chazy-Curzon spacetime are presented 
in Figs. [2]and[3l In Fig. [2] we can see that the shape of 
the effective potential curve depends only on the angular 
momentum £ and that all the curves have two circular 
orbits, except the doted curve. In this graph the doted 
curve has £ — 3.52m and represents the marginally sta- 
ble circular orbit. In the Fig. [3] we show the potential 
with £ — 4.5m and, with the dotted lines, different val- 
ues of the energy E 2 in order to analyze the possible 
orbits. Thus, for E\ we have three different radius as 
the horizontal line cut the potential in three points. For 
the least radius A, the particle is confined in potential 
well as in Fig. [TJ < p/m < A. Whereas for the oth- 
ers radius, C and D, we obtain a bounded orbit between 
these radius, i.e., if we take E\ — 0.975, we obtain a 
bounded orbit between C = p/m sa 7.66 (perihelio) and 
D = p/m w 67.52 (aphelio), Fig. [4] When the energy is 
Ei we get a turning point, so for Ei — 1.02 the turning 
point is B = p/m « 5.45, we show in Fig. [5] the corre- 
sponding graph. For the point F with energy E2 = 1.02, 
the particle is confined in the potential well, as in A. 
Whereas for an energy E 3 , greater then the maximum of 
the potential, there are not turning points and the par- 
ticle moves only in a direction. Finally, in the Fig. [6] we 
present the range of stability for particles moving in a 
circular orbit by plotting the specific angular momentum 
(|3"8")) as a function of the radius of the circular orbit. The 
range of stability is 3.52m < £ < 00 and 5.23m < p < 00. 
The point have coordinates (£/m,p/m) = (3.52,5.23), 
and corresponds to the radius of the marginally stable 
circular orbit. 



IV. SOLUTIONS IN OBLATE SPHEROIDAL 
COORDINATES 



In the oblate spheroidal coordinates, the solution of the 
Laplace equation @ is 



Wn = -^C 2fe P 2fc (7 7 )7 2fe + 1 Q 2fc (iO, 



(42) 



fc=0 



where Cik are constants, Pk are the Legendre polyno- 
mials and Qk are the Legendre functions of second kind 
p^ |. This solution repesents the exterior Newtonian po- 
tential for an infinite family of axially symmetric finite 
thin disks, recently studied by Gonzalez and Reina [2fj| . 
and whose first member, n — 0, is the well-known Kalnajs 
disk [2l[ . We also studied in another paper [22| the kine- 
matics around the first four members of this family by 
means of the Poincare surfaces of section and Lyapunov 
characteristic numbers and fond chaos in the case of disk- 
crossing orbits and completely regular motion in other 
cases. 
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The constants C 2 fc appearing in (l42j) are given by 



mG 
~2a 



rV2 



(4Jfe + l)(2n + l)! 



2 2 ™(2/c + l)(n - k)W(n + k + |)g 2 fc+i(0) 



where <?2fc(£) = * 2/c+1 Q2fe(iC)j 771 i s the mass of the disk 
and G is the gravitational constant. Now, due to the 
presence of the term (n — k)\ at the denominator, all the 
Cik constants vanish for k > n. The variables i] and 
£ are the oblate coordinates related with the cylindrical 
coordinates by 



p 2 =a 2 (l + e)(l-V 2 ), z = a£ri, 



(43) 



where a is a constant, —1 < r) < 1 and < £ < oo. In 
(|43p a is the radius of the thin disk, henceforth we take 
a = 1. In the plane z = 0, we have two regions: if £ = 
then r\ = y/l — p 2 , whereas if 77 = then £ = \J p 2 — 1 . 
These two regions correspond to the regions inside and 
outside of the disk, respectively. 

We now write the different equations for this two re- 
gions in the simple case of null geodesies. So, from (Ti~6"l) . 
the radius of the circular orbits inside the source is 

n 

$>fcC 2fc g 2fc (0) [P2k-i(v) - r]P2k(v)] = V, (44) 
fe=l 



where r\ = y/l — p 2 . The stability condition in oblate 
spheroidal coordinates, inside of the disk, take the form 



77(1 - r] 2 )i>n 



(l+V 2 )^r, < 0, 



that, for the solution (|4"2")l , can be written as 

m 

^2fcg 2fc (0) [2kr]P 2k (r]) +F 2fe _ 1(77)] > 0. 



(45) 



(46) 



fc=i 



Now, for null geodesies outside of the disk, the radius of 
the circular orbit can be obtained from the expression 

n 

2C 2fc P 2fe (0) [£<fefc(0 + <Z2*+i(0] + e = 0, (47) 

fc=0 



where £ = ^/p 2 — 1. The corresponding stability condi- 
tion is 



that using (|42)) . becomes 



(£ 2 - < 0, 



(48) 



J2 C2k(2k + 1) [2£g 2fc (0 + q 2k+ i(0} > 0. (49) 



k=0 



0, was obtained 



The first solution of (|42l) . when n 
independently by Zipoy [23| and Vorhees [24( , and inter- 
preted by Bonnor and Sackfield [25[ as the gravitational 
field of a pressureless static thin disk, this disk is singular 



at the rim. The function ip for the first three members 
of family of disks (g2} is given by [U HI] 



^ = -^[cot-^ + ^^-l)], 

^2 -- 



M[cot- 1 £ + — (3r? 2 -l) 



(50) 



(51) 



V>3 



with 



.4 

c 



+ £(35r7 4 - 3077 2 + 3)], 

-p[cot-^ + ^(3 V 2 -l) 
6 

+ ^(35r/ 4 -30ry 2 + 3) 

+ C(231r] 6 -315?; 4 + 105?7 2 -5)], (52) 



i [(3£ 2 + l)cot- 1 £-3£], 
^[(35£ 4 + 30£ 2 + 3) cot- 1 £ - 35£ 3 - 
[(231£ 6 + 315£ 4 + 105£ 2 + 5) cot" 1 £ 



8448 



231 

-231£ 5 - 238£ 3 - — £], 
5 



wher p = m/a and n = 1, 2, 3, respectively. 
For n = 1, the other metric function 

71 = 9 M 2 (?? 2 - 1) [9£V -£ 2 + 4t7 2 + 4 

2 , i\/nf2 2 t 2 1 2 i\/ ,-1 t\2 



rely. 

is B mil 



+(r + i)(9cv - r + »r - ^M- 1 er 

-2£(9£V - £ 2 + 7r/ 2 + 1) cot" 1 (\ /16. (53) 

This disk is also singular at the rim 28s]. For n = 2, the 
metric function 7 is given by 

7 2 = 25/i 2 (t7 2 - 1) {9£ 6 (1225/7 6 - 1275?7 4 + 3157 ? 2 - 9) 
+3r/ 4 (5050?/ 6 - 3630r? 4 + 366?/ 2 + 6) £ 4 
+ (4945?7 6 - 723 - Abff - 81) £ 2 
-6 [375r? 6 + II377 4 + I577 2 + 27 
+3 (122577 6 - 127577 4 + 31577 2 - 9) <f 
+ (6275?7 6 - 4905?7 4 + 68I/7 2 - 3) £ 4 
+ (3005ry 6 - llllr/ 4 - 105r; 2 + 3) £ 2 ] ^cot" 1 £ 
+9(£ 2 + 1) [%rf + hrf - 5?7 2 - 9 
+£ 6 (122577 6 - 1275r? 4 + 315/7 2 - 9) 
+£ 4 (127577 6 - 785?7 4 + 17?/ 2 + 5) 
+£ 2 (315^ 6 - 1777 4 - 47r! 2 + 5)] (cot" 1 £) 2 
+256 (t] 6 + r] A + r/ 2 + l) } /2048. (54) 

For n > 2, although the metric function 7 can be easily 
obtained by integrating the equations Q - ((4]) properly 
written in oblate spheroidal coordinates, they are not 
explicitely presented here due to their higly involved ex- 
pressions. 

Now we analyze some examples. If we take n = 1, 
the second member of family of disks, we obtain for the 



6 



radius of circular orbit, the angular velocity and specific 
energy corresponding to this radius, the relations 



P 2 = 



3-7T/1 ' 



E 2 
I 2 



2 



(55) 



where /i > 2/3ir, £ is an arbitrary constant and the ra- 
dius corresponds to a stable equilibrium of the effective 
potential for a null geodesic. Now, outside of the source 
the effective potential increases until a maximum value, 
corresponding to the unstable circular orbit, and then 
diminishes until when p increases, according to (|13|) . 
Again, we cannot find the marginally stable circular or- 
bit in this disk for the massless particles. 

The behavior of the effective potential for different val- 
ues of parameters is similar to that we displayed in Fig. 
[51 corresponding to the timelike geodesic. In the effective 
potential of the Fig. [5] we put £ = 8, 6 (gray curve) and 
3.56 (dotted curve). The behavior for different values of 
the parameters is similar, so we take [i = 1 as an exam- 
ple. In the exterior case we consider the third member of 
the family, n = 2, Fig. [8] In this graph the point B cor- 
responds to the radius of the marginally stable circular 
orbit, which have I « 3.688. The points C and D have 
angular momentum I — 4.2 and effective potential 0.935, 
the motion is bounded between the radius 7.61 and 19.14, 
Fig. [9l In this graph we choose 



V. SOLUTIONS IN PROLATE SPHEROIDAL 
COORDINATES 

The general static axisymmetric vacuum solution for 
tp in prolate spheroidal coordinates (it, v) is given by [29| 



^ = ^(-l) n+1 rfnQnWPn(v), 



(56) 



n=0 



where u > 1, — 1 < tj < 1, and the d „, a re constants 
related with the multipole moments 0, [3(| ■ P n are the 
Legendre polynomials and Q n are Legendre functions of 
second kind. These coordinates are related with Weyl's 
canonical coordinates by 



being m the mass of the source that produce the field. 

The asym pto tically flat solution for 7 was found for 
Quevedo in [2{j. The monopolar solution, 1 = 0, with 
do = 1 corresponds to the Schwarzschild field. The solu- 
tion of equation (|16j) for the monopolar solution is u = 2, 
that is the unstable radius of a null circular orbit in the 
Schwarzschild field. For a complete study of motion in 
Schwarzschild field in the equatorial plane, see [3l|. The 
solution of (1561) for I = 2 is the Erez- Rosen metric 1321 , 



Vi 



do u-1 
-7T In 7 

2 u+1 



(58) 



1) 



4^ 



l)ln- 



1 



2 U 



u + 1 

In this solution, g?o and di are related with the monopole 
and arbitrary quadrupole moment, respectively [3(j. The 
study of orbits in this solution was developed by different 
authors 

Now, in this section we expose the expressions for the 
different quantities corresponding to the motion of a par- 
ticle in a circular orbit in prolate spheroidal coordinates. 
That is, the specific energy, the angular velocity and the 
radius of the marginally stable circular orbit, which are 
obtained through of the effective potential in the equa- 
torial plane, 



V = e 2 ^ e 2 + 



(u 2 — 1) 777/ 



p 2 = 777 2 (l-7; 2 )(77 2 -l), 



(57) 



We begin with the equations corresponding to null 
geodesies, when e = 0. So, the radius of the circular 
orbits can be determinate by 

2(u 2 -l)^,„ = «, (59) 
that in terms of (|56p takes the form 
l 

2d 2 nP 2 , i (0)(7 i 2 - l)Q' 2n {u) = u, (60) 

with the stability condition 

u{u 2 - 1)V>,„ U + (u 2 + lty.t, > 0, (61) 
that using (|56[) can be written as 



d2nP2n(0) [(u 2 - l)Q' 2n (u) - 27777,(277 + l)Q 2n (u)] > 0, (62) 

n=0 

I 



where the equality corresponds to the equation for the 
radius of the marginally stable circular orbit. Finally, 
the other expressions are found by means of (jT4"]) and 



On the other hand, for the massive particle we obatin 
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the expressions 



E 2 = 



e~ 2 ^ (u 2 -l) 2 m 2 ^ u 
u-2( U 2 -l)^„ ' 

e 2V [ u _ (u 2 _ 

u-2(u 2 -!)</>,„ ' 



(63) 
(64) 



where u — 2 (it 2 — l) V\« > in order that the energy per 
mass unit and the angular moment per mass both be not 
imaginary. The stability condition and the radius of the 
marginally stable circular orbit are 

VV {3u 2 + 2 + 2(u 2 - 1)V>,„ [2(u 2 - - 3u] } 

+ u{u 2 - l)ij) ;UU > 0, (65) 

with u - 2 (u 2 - l) f/v > 0. 

Finally, we show in Fig. ITOthe region of stability by 
means of a graph of the specific angular moment (|63[) . In 
particular, by taking arbitrary values of the parameters 
in the Erez-Rosen solution, one can see that the range is 

< u < 4.77, 4.77 < u < oc, 

and 11.62 < £ < oo for stable orbits, where we choose 
da = 1 and di = 1.5. Here the marginally stable circular 
orbit have coordinates (£/m,u) ~ (11.62,4.77). 



VI. CONCLUDING REMARKS 



all of them obtained through of an effective potential. 
The specific angular momentum was presented for the 
timelike geodesic and was used to determine the range of 
stability of the orbit of the particle, so the the minimum 
value represents the marginally stable circular motion. 

In order to find the trajectory of the particle, we ana- 
lyzed the analytical results obtained. The character of 
the motion is determined essentially by means of the 
behavior of the effective potential. Thus, we displayed 
different graphs of effective potential before resolve the 
differential equations of motion of the particle. Then, 
we began with the Chazy-Curson field in the case of the 
cylindrical coordinates, as discussed in section III. The 
motion of particles around oblate deformed bodies was 
developed in section IV, by means of the analysis of the 
properties of some member of the family disks (1421) . On 
the other hand, the prolate case was presented in section 
V, where we found the range of stability of the Erez- 
Rosen solution in the special case of mass particles. 

In summary, we concluded that for these solutions all 
the circular orbits are unstable only in the case of null 
geodesic, whereas that do not exist marginally stable cir- 
cular orbit for null geodesic. In contrast, we found that 
for mass particles the orbits can be unbounded, bounded 
or circular. This behavior can be seen by means of the 
effective potential and verified through the numerical so- 
lution of the equations of motion. Moreover, for the time- 
like geodesic we found the radius of the marginally stable 
circular orbit in different coordinate systems, cylindrical, 
prolate and oblate. 



In this paper we analyzed the behavior of free test 
particles in the equatorial plane of static axisymmetric 
spacetimes. We presented several general expressions for 
the circular orbit in null and timelike geodesies: radius, 
specific energy, specific angular momentum, angular ve- 
locity and radius of the marginally stable circular orbit, 
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FIG. 2: Effective potential as a function of p/m for timelike geodesies, with different values of i/m, for the Chazy-Curzon field. 
For the curves from top to bottom we take i/m = 4.1,3.8,3.65 and t = 3.52m, respectively. The point in the dotted curve 
corresponds to p/m = 5.23, which is the radius of the marginally stable orbit with i/m = 3.52. 
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FIG. 3: Effective potential for timelike geodesic with I — 4.5m in the Chazy- Curzon spacetimes. 




FIG. 4: Orbit of the particle corresponding to the effective potential of the Fig. [3] with E\ = 0.975. This orbit is bounded 
between C w 7.66 and D w 67.52. The initial conditions are p(t = 0) w 0.17, ip(t = 0) = tt/6 and p(t = 0) = 20. 
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FIG. 5: Orbit of the particle corresponding to the effective potential of the Fig. [3] with E2 = 1.02. The initial conditions are 
p(t = 0) « 0.17, if = tt/6 and p = 20. 
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FIG. 6: The specific angular momentum, i/m, as a function of the radius of the circular orbit, p/m, for timelike geodesic in 
the Chazy-Curzon field. The point have coordinates (£/m, p/m) = (3.52, 5.23). The range of stability is 3.52m < £ < 00 and 
5.23m < p < 00. 
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FIG. 10: Specific angular momentum, £/m, as a function of u for the circular orbit in the Erez-Rosen field with do = 1 and 
d 2 = 1.5. 



